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Abstract 

In this review, we present our recent computational work on carbon-based 
nanostructured composites. These materials consist of carbon crystallites em- 
bedded in an amorphous carbon matrix and are modeled here through classical 
and semi-empirical quantum-mechanical simulations. We investigate the ener- 
getics, mechano-elastic, and optoelectronic properties of these materials. Once 
the stability of the composites is discussed, we move on to the calculation of 
their elastic moduli and constants, their anisotropy and elastic recovery. At a 
next step, we focus on diamond composites, which were found to be the most 
stable among the composites studied, and went beyond the elastic regime to in- 
vestigate their ideal fracture. Finally, for these materials, the electronic density 
of states, dielectric function, and optical response were calculated and linked to 
the disorder in the structures. Our findings unveil the high potential of these 
materials in nanotechnological applications, especially as ultlra-hard coatings. 
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1. Introduction 

Carbon-based nanostructured composites (na-C) [l|, are believed to in- 
termingle properties of carbon nanostructures with those of single-phase amor- 
phous carbon (a-C)Q. These nanocomposites are hybrid forms of carbon in 
which nanocrystallites are embedded in the a-C matrix. In this sense, their 
properties can be tailored by controlling the type and size of the embedded 
structures, as well as the matrix material. Such composites, also investigated in 
this work, are relevant to both bulk and thin film materials. 

A variety of carbon nanostructures can compose na-C, ranging from diamond 
crystallites, to entirely three-dimensional sp 2 covalent conformations with no 
dispersion/ van der Waals (vdW) bonding, i.e. open graphene structures with 
negative curvature (schwarzites) 0, [H, @) , to carbon nanotubes (CNT) 0| bound 
by vdW forces. Na-C films have been synthesized before by using different 
methods The synthesis of nano diamond composites (nD/q-C) and its growth 
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mechanisms have been explored in both hydrogenated 0, [l3| and pure 
tetrahedral amorphous carbon (ta-C) matrices. 

Here, we review our recent theoretical and simulational studies of nanos- 
tructured carbon films, which shed light on their mechanical and optoelectronic 
properties. This review is organized as follows: Section [2] outlines the method- 
ology followed here. Section [3] discusses the energetics and stability of all mod- 
eled nanocomposites. The elastic and mechanical properties of diamond and 
nanotube composites are investigated in section 0J The ideal fracture and op- 
toelectronic properties of diamond nanocomposites are discussed in sections [5J 
[5J We summarize in section [7] 



2. Methodology 



We have recently followed a theoretical approach to investigate na-C and 
unravel the fundamental principles governing the interaction of nanostructures 
with the a-C matrix (see for example 13|,llJ]). Two different methods were used, 
a Monte Carlo (MC) method within the Tersoff empirical potential approach 
15| . and a tight-binding molecular dynamics (TBMD) scheme which utilizes 



separately two hamiltonians, the Mehl & Papaconstantopoulos [16|, [17| (TB- 
NRL), and the environment-dependent (EDTB) Wang & Ho model. The 
former employs non-orthogonal atomic orbitals, while the latter goes beyond 
the traditional two-center approximation and allows the tight-binding parame- 
ters to change according to the bonding environment. Using both the MC and 
TBMD, one is able to combine higher statistical precision with higher accuracy 
and bridge the gap between classical and first-principles calculations. Both the 
Tersoff potential, as well as the two hamiltonians have been used to successfully 
model the interatomic interactions in various crystalline and amorphous carbon 
forms [II [11 [II il, [23 1 

The TBMD simulations are carried out in the canonical (N, V, T) ensemble. 
The volume V and the number of particles (atoms) N are constant. Disper- 
sion/van der Waals (vdW) forces within the Tersoff potential approach, will 
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be included by adding a simple Lennard- Jones potential, which successfully de 



scribes properties of carbon nanocrystals [2J, |25| ■ The parameters used for this 
potential are e = 2.964 meV and a = 3.407 A. 

For the calculation of the optoelectronic properties the EDTB scheme was 



used, which has proved to be quite efficient in similar cases [22fl . The electronic 
structure can be easily extracted and is directly correlated to the imaginary 
part £2 of the complex dielectric function. This is in turn yields the optical gaps 
and absorption coefficient, which further link to the disorder in the structures 
through the Urbach energy. Additional details can be found elsewhere (26j . 

2.1. Structure generation 

We have generated the nanocomposites studied here by melting and sub- 
sequent quenching at constant volume various carbon crystalline cells, while 
keeping a certain number of atoms in the central portion of the cells frozen in 
their ideal crystal positions. Consequently, the cells are thoroughly relaxed with 
respect to atomic positions and density. For the nanotube composites, a core 
is removed from a diamond crystal leaving a void with a shape similar to that 
of the nanotube, in which the nanotube is then placed. Details can be found in 

Ref. dEG|. 



The total number of atoms in the structures ranges from ~ 2500 to 4536, 
and the number of atoms in the nanocrystals ranges from ~ 300 to 900. The 
size of the computational cells ranges from about 1.5 to 4.0 nm. For the TBMD 
relaxed structures we use smaller cells up to 512 atoms. The diameter (d) of all 
crystalline inclusions used are in the range d=0.6-2.8 nm. Periodic boundary 
conditions (PBC) are applied to the cells. The exact details for each composite 
studied with the Tersoff potential are summarized in Table [TJ Diamond inclu- 
sions are almost spherical, while CNTs have open ends and extend through the 
entire length of the cube. Due to the PBC, this corresponds to a dense array 
of CNTs packed in parallel, which is an idealized model of a CNT nanocom- 
posite. However, this model captures all essential features of the CNT-matrix 
interaction resembles CNT bundles with a-C material in between the tubes. 

Representative structures of a specific group of conformations are chosen and 
embedded in amorphous matrices: the sp 3 conformations (a) Diamond (D), and 
the (b) high-pressure phase BC8 structure for carbon, the negatively-curved sp 2 
conformations with a D periodic minimal surface, (c) polybenzene (6.8 2 D) [J], 
(d) PCCM (porous conducting carbon modification), an sp 2 -bonded structure 
with eightfold rings [27[, (e) C168, a low-density, negative-curvature schwarzitc 
structure [Bj], (f) the positively curved sp 2 carbon nanotubes (CNT) Q with 
single (SWNT) and multiple (MWNT) graphitic walls. Structures (c)-(e) do 
not involve vdW bonding. In FigfTJ representative nanocomposites of the ones 
studied here are shown: a diamond of 1.70 nm diameter, a C168, and a (7,10) 
SWNT of 1.17 nm diameter, embedded in a highly tetrahedral (z am ~ 3.8, 
p ~ 3.3 g cm' 3 ) and low density (z am ~ 2.8, p ~ 1.5 g cm' 3 , and z am ~ 
3.2, p ~ 2.2 g cm' 3 , respectively) matrices. The mean coordination of the 
amorphous matrix and the total density of the composites are denoted z am and 
p, respectively. FigQJb) shows clearly, that no covalent bonding between CNT 
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atoms and the matrix exists, as the vdW forces drive a reconstruction of the 
surrounding medium, producing a graphitic wall at a distance of 3.4 A from the 
CNT. 

3. Formation Energies 

We begin with the investigation of the stability of the nanocomposites stud- 
ied here. The interaction of the embedded configuration with the host is taken 
care of by defining the formation energy of a nanocrystal, given by 

Eform = Efotal ~ N a E a — N C E C (1) 

In this equation, E to tai, E c are the total cohesive energies per atom of the 
composite system and the crystalline phase, respectively. N c , N a are the number 
of atoms in the nanocrystal, and the amorphous matrix, respectively, and E a 
is the cohesive energy per atom of the pure amorphous phase (without the 



nanocrystal) with coordination z am l_3j. A negative value of Ef. 



stability of the nanostructure, a positive value indicates metastability. 

The formation energies, hence stabilities, of all nanocomposites in this study 
are depicted in FigJ^l Of the nanostructures considered, only diamond and 
the Ci68 schwarzite are stable (Ef orm < 0), the former at high z am (p am > 
3 <?cm -3 ), the latter at low z am (p am < 1-2 gcm~ 3 ). Both crystallites are 
stable in a matrix with structural characteristics similar to their own, which 
relates to a small density gradient between the matrix and the inclusion: dia- 
monds(schwarzites) are dense sp 3 (low density sp 2 ) networks, thus are embedded 
in a dense(dilute) matrix. A small density gradient does not give rise to dis- 
continuities and thereby energetically favors the embedding of the crystalline 
inclusion. Overall, diamonds are more stable than schwarzites. Their enhanced 
stability in highly tetrahedral networks indicates that once nucleation centers 
are formed large diamond regions can be further grown, as shown in experiments 
@, EH- Through annealing of the diamond composites, an unstable structure 
extensively shrinks, and only a small core remains intact, while the stable struc- 
ture retains its shape and should expand. Metastable diamonds are on the 
average larger, but more deformed, than stable diamonds in ta-C. For a better 
nanostructured material it is preferable to nucleate diamond cores in a ta-C 
matrix. 

All other nanostructures studied in this work are metastable (Ef orm > 0) 
through the whole coordination range. Their curves exhibit a well defined local 
minimum at which Ef orm is quite small, suggesting that they can be main- 
tained in the amorphous matrix under moderate conditions and their synthesis 
is possible. The fact that Ef orm for 6.8 2 D and PCCM composites are only 
slightly lower than that for a diamond inclusion suggests the possibility of a 
phase transformation of sp 2 -bonded structures to diamond under the proper 
conditions of pressure and temperature. This would subsequently aid the ex- 
periments towards synthesis of such theoretical carbon crystals. More details 
are discussed in [l3| . CNT composites are also metastable with minima in their 
Eform at intermediate coordinations z a 
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4. Elastic and mechanical properties 

The elastic and mechanical properties of carbon composites are important 
for potential applications as ultra-hard protective or biomedical coatings. We 
focus on diamond nanocomposites (nD/o-C), which are the most stable among 
the ones we dealt with in this work, and the nanotube nanocomposites (CNT/a- 
C), as a representative sp 2 structure with positive curvature which is of wide 
technological interest. 



4-1- Bulk and Young's moduli 

We begin with the calculations of the bulk modulus (B), which was calcu- 
lated by fitting the energy vs volume curve after a hydrostatic deformation to 
the Birch- Murnacham equation [28[. Figs. EJa) , HJ^a) show the calculated mod- 
uli plotted as a function of the nanocomposite density for various diamond and 
nanotube diameters, respectively. The bulk modulus is significantly enhanced 
with increasing density and inclusion size. Replacement of some amorphous ma- 
terial by crystalline increases noticeably the bulk modulus, with respect to B of 
a-C, in both diamond [23[ and nanotube composites (FigfJJa)). B for some sam- 
ples is even higher than that of the amorphous diamond network (WWW) [2^ ] 
and close to that of diamond 23]. A similar enhancement of the hardness with 
the grain size was found in ultra-nanocrystalline diamond [3Cj | and copper films 



3l| . A comparison of B for S WNT and MWNT inclusions of the same diame- 



ters embedded in a-C of similar densities reveals, that the number of walls that 
consist a MWNT is crucial for the whole composite's rigidity (see FigfJJb)). By 
increasing the number of graphitic walls within a CNT the composite becomes 
more rigid. 

A check whether the rule of mixtures [32[ applies to diamond and CNT 
composites revealed that it indeed holds for the former, but not for the latter 
structures. Apparently, the CNT composites are significantly influenced by the 
presence of the matrix. In these structures, the microstructural details of the 
matrix have been significantly altered with respect to pure a-C. We note again, 
that the effect of the inclusion is to drive the surrounding matrix atoms further 
away from the nanotube and stabilize them at a distance of about 3.4 A. On 
the other hand, no such a significant effect was evident in nD/a-C, as there is a 
smooth structural transition from the inclusion to the matrix. In other words, 
there is a direct connection of the inclusion-matrix bonding and the validity 
of the rule of mixtures: in the nanotube composites the bonding between the 
matrix and the inclusion is of a van der Waals type and the rule does not hold, 
while in the diamond composites the two phases are covalently bonded and this 
rule still holds. 

The trends for the Young's moduli (Y), which is easily extracted from the 
elastic constants discussed next, of both composites are the same as in the case 
of B. For nD/a-C, Y is shown in FigJ3] A considerable enhancement of Y 
was found, especially in high density matrices, with values approaching that 
for diamond [23j. In the case of CNT/a-C, a marked characteristic is the high 
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elastic anisotropy which is monitored by the ratio 

A = Y ± /Y {1 (2) 

where Yii , Yj_ are the components of Y in the direction of the tube (axial) and 
in the transverse directions. A composite with a SWNT with diameter of 1.65 
inn and a density of 1.99 g cm~ 3 (2.14 g cm~ 3 ) has Y\\ = 478(564) GPa and 
Y± = 231(306) GPa, yielding A = 0.48(0.54). The anisotropy increases with 
increasing tube diameter, as the tube contribution overwhelms the isotropic 
matrix part. 



4-2. Elastic constants 

To supplement the calculation of B, the remaining elastic moduli are com- 
puted. All simulation cells of the diamond composites are cubic systems, thus 
only the knowledge of the en, C12 and C44 is needed. This is not exactly the case 
for CNT composites, as these are modeled through tetragonal computational 
cells. To calculate the elastic constants, we apply the appropriate deformation 
to the system and compute its total energy as a function of the imposed strain. 
The curvature of this function at its minimum yields the desired modulus. The 
results are plotted in Fig(3Jb) as a function of the density for diamond compos- 
ites. All elastic constants increase with increasing density. This is also the case 
in CNT/a-C, for which representative values are summarized in Tabled For 
an isotropic material, 

?? = en - C12 = 1 ^ 

2C44 

which holds within 4% or less, for nD/a-C. For a typical CNT composite, 77 = 
0.75, while for all pure a-C networks i] — 0.98. 

We have also looked at the elastic recovery during a compression-decompression 
cycle from to 100 GPa for diamond and nanotube composites (with inclusions 
of sizes 1.7 nm and 1.2 nm, respectively), as well as for amorphous structures 
[33] ]. The elastic recovery of nD/a-C and single-phase a-C is almost perfect 
in the highly tetrahedral regime (ta-C), but is not perfect as the sp 3 fraction 
declines. An almost fully sp 2 -bonded a-C network regains its initial volume at 
about 95 %. A CNT composite, having a high sp 2 fraction does not fully recover 
and shows a hysterisis in the compression curves. Our results indicate that the 
compressibility from CNT/a-C towards a-C and nD/a-C is decreasing, as long 
as the sp 2 fraction remains high. 



5. Ideal strength of diamond nanocomposites 

A study beyond the elastic regime gives an insight to the ideal fracture of 
a-C and nD/a-C Tensile (along diamond's easy slip {111} plane) and shear (on 
the {111} plane in a (112) direction) loads were applied and the stress versus 
strain curves were extracted. The NRL hamiltonian at low temperature was 
used in this study. This method has already been tested in the case of diamond 



and the findings compared well with previous ab initio results [34L |3 
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Three a-C samples were considered (a typical non-hydrogenated a-C with 
z am —3-4:7, a ta-C sample with z om =3.78 and a WWW sample) and a nD/a-C 
with an inclusion of about 1.2 nm embedded in ta-C with an sp 2 fraction of 
about 30%. The stress vs. strain curves in the case of the tensile load are 
summarized in FigJSJ left panel. The maximum stresses are roughly 60, 40, 30, 
and 124 GPa for z am =4.00, 3.78, 3.47, and diamond, respectively. Thus, the 
ideal strength of a-C is proportional to its sp 3 ratio. The stress vs. strain curves 
for the low-density a-C (sp 3 fraction ~ 50 %) sample suggest a ductile behavior. 
The fracture of a-C, which has an ideal strength half that of ta-C, takes place 
at a higher load than ta-C, i.e. larger strain needs to be applied to break the 
a-C network. Similar are the results when a shear load is applied. 

For the diamond composite, the easy slip plane of the embedded diamond 
core was not found to play a crucial role in the total strength and no weak di- 
rections were present. The crystalline phase remains unaffected by the external 
load, which is almost completely taken by the surrounding amorphous matrix. 
The response, thus, of nD/a-C to external load beyond the elastic regime is 
identical to the response of the embedding matrix. The onset of fracture takes 



place at the weakly bonded sp 3 sites in the amorphous matrix 23j. Bonds in 
the crystal are at this point actually elongated. Fracture in diamond nanocom- 
posites is inter-grain and their ideal strength is similar to the pure amorphous 
phase. A snapshot of a nD/a-C at the onset of fracture (i.e. at the maximum 
stress) is shown in FigO right panel. 



6. Optoelectronic properties 

Finally, focus is given on the electronic and optical properties of diamond 
nanocomposites. We were able to probe at a local atomic level the optoelectronic 
response of the composite and link it to the associated disorder, which exists at 
the interface between the diamond inclusions and the embedding a-C matrix. 
This can be quantitatively probed by extracting the Urbach energies from the 
optical parameters. Disorder in the nanocrystals appears in their outer shell 
near the interface and is manifested as bond length and angle distortions. 

As representative quantities of an analysis done elsewhere [26| , the electronic 
density of states (EDOS) and the imaginary part (£2) of the dielectric function 
for nD/a-C are shown in Fig(5J The same approach for the study of the di- 
electric response for diamond- like amorphous carbon has been used before (3r3 | . 
The results compare qualitatively well with recent experiments [37| , taking into 
account that in the latter case the diamond-like film was hydrogenated. A first 
look at the EDOS in Figj6] reveals that the main contribution to the forma- 
tion of the gap, i.e. the absence of electronic states, and of any states located 
in the gap region comes from the amorphous matrix part, as this has a lower 
energy gap than diamond due to the presence of sp hybrids. As the sp 2 ra- 
tio, which controls the gap through the separation of it-it* bonding-antibonding 
states increases, the gap region shrinks. A smaller contribution to the gap states 
originates from distorted sp 3 atoms of the diamond inclusion at the interface. 
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Decomposition of the total £2 into contributions from the atoms in the crys- 
talline and amorphous regions showed, that the optical response of the nanocom- 
posite system is progressively varied. This variation occurs as the sp 2 fraction 
in the amorphous part is increased from having nearly diamond-like nature to 
having nearly graphite-like nature. It was also clear that the energy gaps, i.e. 
absence of electronic transitions at the edge of the spectra, are reduced as the sp 2 
fraction increases. Near the gap the main contribution comes from transitions 
between states localized in atoms of the amorphous part. 

Calculation of the Urbach energy (Ejj) in nD/a-C, which is linked directly 
to the optical response provides information about the disorder in the structure. 
We have recently shown [38] . that is a good measure of disorder in a-C, as 
it reflects not only bond-length and bond-angle distortions (structural disor- 
der), but also topological disorder related to the size distribution of chains and 
clusters of sp 2 and sp 3 atoms in the mixed phase Thus, it is a measure of 
the inhomogeneous disorder in the two-phase network. For example, in highly 
tetrahedral ta-C, the larger the sp 2 cluster sizes the higher the Ejj values. This 
is because the embedding of minority configurations into a host phase produces 
disorder, which is reflected in the tt-tt* subgap transitions. We also showed [Hf 
that transitions between states localized on unpaired sp 2 sites (dangling bonds), 
or between such states and neighboring ir or a states, are negligible. This is in 
contrast to a-Si:H, where dangling bond defects are the main contributors to 
the Urbach edge. 

E[/ for a typical composite is shown in Fig|7]and underlines a nonmonotonic 
behavior for both the composite and pure a-C. The variation of Ejj in the case 
of nD/a-C comes from the atoms in the matrix, but with higher values at the 
maximum of the variation and at the tetrahedral region, indicating the excess 
disorder produced on these atoms. We have also broken Ey down into atomic 
components by calculating separately the absorption coefficient (a) contributed 
by the atoms close to the diamond nucleus and close to the interface. Results 
shown in Fig indicate that the variation in a nearly coincides with the contri- 
bution from the atoms at the interface, while the absorption from the nucleus 
approaches the absorption of crystalline diamond. The atoms at the interface 
determine the Ejj for the whole composite. Defects at the interface (unpaired 
sp 2 sites and dangling bonds) contribute negligibly to the optical transitions and 
the Urbach edge, as in pure a-C 38]. On the other hand, Ejj for the diamond 
nanoinclusions takes quite high values (~ 0.6) and it is nearly constant as a 
function of sp 3 fraction. This shows the significant structural disorder at the 
outer shell of the nanocrystals. 

7. Summary 

In this review, we have briefly presented some of our work on the energetics, 
stability, mechanical and optoelectronic properties of carbon-based nanostruc- 
tured composites. These are materials consisting of a carbon crystalline nanoin- 
clusion embedded in an amorphous carbon matrix. The inclusion ranges from a 
diamond crystallite to sp 3 and sp 2 crystalline conformations with negative and 
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positive curvature. The ability to control the mechanical and optoelectronic 
properties of these materials by varying the density of the embedding medium 
and the type and size of the crystalline inclusion opens up enormous possibilities 
for potential technological applications of carbon nanocomposites. 

We found diamonds and schwarzites to be highly stable in dense amorphous 
carbon matrices, while the other sp 3 and sp 2 conformations are metastable. A 
considerable enhancement of the elastic moduli for nD/a-C and CNT/a-C was 
found as compared to pure a-C values with increasing density and inclusion size. 
a-C and nD/a-C of high density showed perfect elastic recovery, which was not 
the case for CNT/a-C. Beyond the elastic regime, nD/a-C fracture inter-grain 
and in a manner similar to pure a-C. Regarding the optoelectronic properties 
of nD/a-C, the matrix was shown to have a dominant role in forming both 
the EDOS and the dielectric function, with increasing sp 3 fraction and density. 
Analysis of the optical response and the Urbach energy, showed that much of 
the disorder in the composites originates from the interface region. 
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Tabic 1: Computational details of the composites with the various inclusions presented in this 
review. Ntot, are the total number of atoms. N cr , and N am are the number of atoms of the 
crystalline inclusion and the amorphous matrix, respectively, d and A arc the initial inclusion 
diameter and the filling ratio. Note, that not all inclusions are spherical; for these, d gives a 
rough estimate of their size. Next to the CNTs the chirality indices (n,m) are given. 



Table 2: Elastic moduli (in GPa) for SWNT nanocompositcs of different densities p (in gem 3 ) 
and nanotubc diameters d (in nm). 
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List of Figure Captions 



Figure 1: Ball and stick models of carbon nanocompositcs: (a) diamond, (b) schwarzite Ci68> 
and (c) a (7,10) SWNT, embedded in a-C matrices of high (for diamond) and low (for Cies 
and the SWNT) densities. The atoms of the matrix and the embedded crystallite arc indicated 
by white and black circles, respectively. 



Figure 2: Formation energies for all carbon composites studied here. Panel (a) shows the 
na-C, apart from the SWNTs which are shown in (b) for various diameters denoted in the 
legend. Labels are self-explaining according to the text. Lines are fits to the points. 



Figure 3: (a) Bulk (B) and Young's (Y) moduli, and (b) elastic constants for nD/a-C as 
functions of the total density (p). The average diameter of the diamond inclusion is 1.7 nm. 



Figure 4: (a) Bulk moduli for CNT/a-C of various diameters (shown in the legend) as functions 
of the density of the composite. The dotted line is a fit to the a-C data, (b) Comparison of 
B for SWNT and MWNT composites with diameters of 2.06 nm. 



Figure 5: Left panel: Stress vs strain curves for structures under tensile load along the {111} 
direction (see text for details). Right panel: Ball and stick model for fracture in the nD/a-C 
for which the stress vs. strain curves are shown on the left. Atoms are shaded according to 
the legend. Large spheres (with label (b)) denote atoms with broken bonds. The snapshot 
corresponds to the maximum stress, at which bonds have started to break. 



Figure 6: EDOS (left) and calculated E2 (right) of nD/a-C with sp 3 fractions and densities (a) 
88 % and 3.24 g cm~ 3 , (b) 71% and 2.91 g cm" 3 , and (c) 51% and 2.58 g cm~ 3 , respectively. 
The black, blue, and red lines refer to the total function, and the inclusion and a-C matrix 
projected contributions, respectively. The Fermi level is at eV. 



Figure 7: (a) Variation of Ejj in nD/a-C as a function of the sp 3 fraction in the matrix (AM). 
Filled circles and squares are the contributions from the inclusion and the matrix, respectively. 
Open circles are values from core atoms in the inclusion. Diamond and triangle are values for 
diamond and amorphous diamond, respectively. Light gray line shows Ejj in pure a-C. (b) 
Breakdown of a (dashed-dotted red line) of a nD/a-C with 88% sp 3 fraction and 3.24 g cm~ 3 
in the a-C matrix into contributions from the core (black solid line) and the interface (black 
dashed line) atoms. Solid brown line indicates the calculated absorption of diamond. 
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